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1.0  Introduction 

In  many  digital  signal  processing  systems,  it  is  necessary  to 
convert  an  array  of  quadrature  saiqples  to  an  envelope  array,  a  process 
equivalent  to  rectifying  in  analog  systems.  For  eetch  quadrature  sample, 
the  quantity  to  be  ccnputed  is 


where  x  and  y  are  the  in  phase  and  quadrature  components  of  the 
6  8 

signal. 

In  practical  systems,  and  y  are  carried  as  Integers  and  the 
output  must  also  be  an  Integer.  Also,  It  is  desirable  to  avoid 
overflow  problems  associated  with  squaring,  and  performing  a  square 
root  is  often  prohibitive  due  to  computation  time-data  rate  treuleoffs. 
A  slnple  estimate  of  r  which  avoids  squaring  and  square  roots,  but 
which  does  not  cause  errors  that  would  affect  the  performance  of  the 
system,  is  therefore  required. 

This  report  describes  this  goad  by  using  linear  combinations  of 


X  *  max(|x^|),|y^|),  y  =  mindx^l,  |y^|  ) 


(1-2) 


of  the  form 


9  »  ax  +  by 


(1-3) 


to  compute  an  approximation  or  estimate  of  r  as  defined  by  Eq.  (l-l). 
Considerable  attention  is  devoted  to  the  implementation  by  simple 
hardware  or  software  oriented  algorithms. 

The  evaluation  of  approximations  of  the  type  in  Eq.  (l-J)  offers 

many  conceptual  difficulties.  In  order  to  minimize  the  complexity  of 

the  concepts  and  mathematics  required  to  analyze  the  approximations, 

they  are  presented  in  terms  of  a  new  set  of  variables.  For  exan^le, 

in  the  (x  ,y  )  or  (x,y)  plane,  r  as  given  by  Eq.  (l-l)  is  no  single 
s  s 

point  or  curve  but  is  a  set  of  circles  centered  at  the  origin.  In  the 
(*  /^jy  /r)  or  (x/r,y/r)  plane,  r  becomes  a  single  circle  of  radius  one 
and  centered  at  the  origin,  called  the  unit  circle.  However,  in  this 
plane,  linear  combinations  of  x  and  y,  such  as  f  as  given  by  Eq. 
are  circles  passing  through  the  origin.  However,  the  con5)lex  variable  w, 
where 


V  =  =  (x/f,y/r)  (l-*+) 

and  r  is  a  linear  combination  of  x  and  y  as  given  by  Eq,  (l-5)>  offers 
the  following  advantages:  (a)  in  the  w  plane,  where  curves  are  plotted 
as  Re(w)  versus  Im{w),  r  as  given  by  Eq.  (1-3)  is  a  straight  line,  and 
(b)r  as  given  by  Eq.  (l-l)  is  the  unit  circle  (a  circle  of  radius  one 
centered  at  the  origin)  as  in  the  (x/r,y/r)  plane.  In  the  w  plane,  it 
becomes  obvious  that  the  optimal  estimates  of  r  are  polygons  approxi¬ 
mating  the  unit  circle,  and  other  concepts  are  equally  simplified. 

The  most  important  of  these  concepts  are  (a)  the  use  of  (x,y)  as  given 
by  Eq.  (1-2)  allows  consideration  of  only  half  the  first  quadrant 
rather  than  the  entire  unit  circle,  (b)  the  use  of  n  straight  lines 
in  this  region  results  in  an  8n  sided  polygon  approximating  the  unit 
circle  in  the  w  plane,  and  (c)  the  use  of  r  =  x  (i.e.,  b=0  in  (Bq.  1-5)) 
for  small  y  results  in  an  8n-4  sided  polygon  in  the  w  plane. 
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The  notations  used  throughout  this  report  are  as  follows. 


Quantity 


y 

r 


e 

s 

b 


Definition 

In  phase  sample. 

Quadrature  sample. 

Largest  of  absolute  values  of  x  ,y  . 

6  S 

Smallest  of  absolute  values  of  x  ,y  . 

8  8 

True  value  of  envelope. 

Estimate  of  envelope. 

Peak  absolute  error  of  estimate, 
rms  error  of  estimate. 

Multiplicative  bias  of  estimate. 
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2.0  Basis  of  Solutions 


i.e.,  the  signal  phase  6  is  uniformly  distributed  over  [-«, n], 

s 

If  the  absolute  value  of  both  samples  is  taken  and  x  is  arbitrarily 
taken  to  be  the  largest,  then 


P(0) 


0  5  e  s  ^ 

otherwise 


(2.4) 


i.e.,  the  modified  signal  phase  6  is  uniformly  distributed  over  [O, jt/4]. 


2.2  Coordinate  Systems 

The  most  obvious  coordinate  system  for  consideration  of  the 

problem  at  hand,  namely,  the  x  ,y  plane,  has  the  fundamental  drawback 

s  s 

that  in  this  system  r  is  a  family  of  circles  rather  than  a  single 

locus.  Therefore,  the  x  >y  plane  coordinate  system  is  not  used  in 

s  s 

this  report. 

The  y/x, r/x  plane  constrains  the  envelope  to  the  locus 

r/x  =  yjl+iy/xr  ,  0  5  y/x  s  1  , 

but  the  weakness  of  this  system  is  that  y/x  has  the  probability 
density 

p(y/x)  =  g  ,  0  S  y/x  S  1 

i+(y/x)‘^ 


a  form  which  ccm^licates  evaluation  of  b  and  s.  Also,  this  coordinate 
system  is  related  to  optimal  estimates  by  conqplicated  relationships. 


6 


The  x/r,y/r  plane  constrains  the  envelope  value  to  the  unit 
circle.  Estimates  of  the  form  f=ax+by  in  normalized  form. 


r/r  =  a  cos©  +  b  sinfl 


,  cosS  =  x/r  ,  sinS  =  y/r  ,  (2-5) 

/  2  2\1-/2 

are  circles  passing  through  the  origin  with  diameter  (a  -fb  j  and 
centered  at  (a/2,b/2),  as  shown  in  Fig.  1.  The  angle  of  maximum  r,  ♦,  is 


=  tan“\b/a),  f  at  maximum 


(2-6) 


The  Q  of  Eq.  (2-5)  Is  the  same  as  that  of  Eq.  (2-U).  This  is  the 
coordinate  system  in  which  the  solution  will  be  presented. 


Another  coordinate  system  of  lTQ>ortsnce  is  the  w  plane,  where 
w  is  the  coo?)lex  variable 

w  =  r/5  (cosfl+j  sine)  =  u+Jv  =  +  j  .  (2-7) 

The  variable  w  is  the  complex  conjugate  of  the  reciprocal  of  r/r 
expressed  as  a  complex  veuriable. 


1  _  jr _ 

w  "  x+Jy 


I  e3Q>(-je) 


The  true  envel<^  is  agsdn  the  unit  circle,  but  Eq.  (2-5)  in  the 
w  plane  is  a  straight  line, 


au  +  bv  ■  1  ,  (2-8) 

with  slope  -a/b  (or,  tan(^  +  it/2),  with  closest  approach  to  the  origin 

/  2  2  \“  1/2 

I a  +b  )  .  The  line  corresponding  to  the  circle  in  Fig.  1  is  shown 

in  Fig.  2.  The  angle  t  is  the  angle  of  minimum  |w| . 
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3.0  Optimal  Estimates 


Optimal  estimates.  In  the  context  of  this  report,  are  taken  to 

mean  the  use  of  n  stred^t  lines  in  the  w  plane  (l.e.,  simple  linear 

combinations  of  x  and  y  )  to  approximate  the  unit  circle  with  minimum 
8  8 

rms  error  s^.  The  multiplicative  bias  b^^  is  irrelevant  to  signal 
processing  considerations.  In  the  v  plane, 

w  =  u  +  3v  ,  v/u  =  y/x  =  twoB  , 

and  6  is  uniformly  distributed,  so  an  optimal  fit  for  a  fixed  number 

of  lines  will  be  a  regular  polygon.  Independent  of  the  loss  or  error 

criterion.  The  use  of  the  absolute  values  x  and  y,  x  2  y  k  0,  Instead 

of  X  and  y  restricts  9  to  the  range  0  to  n/4,  and  n  llne^  combina- 

tlons  of  X  and  y  are  equivalent.  In  general,  to  8n  linear  combinations 

of  X  and  y  .  If  the  line  at  9=0  In  the  v  plane  is  vertical  {9*cx), 

8  8 

or  if  the  line  at  0=x/h  in  the  v  plane  has  a  slope  of  -l(r*c(xfy)), 
the  lines  In  the  w  plane  will  be  parallel  to  tangents  to  the  unit 
circle  at  those  points  and  the  polygon  will  have  8n-4  or  8h-8  sides. 
Thru,  the  ^c(x4-y)  case  Is  to  be  avoided;  however,  the  simplicity  of 
^cx  near  9>*0  makes  this  case  worthy  of  consideration.  Accordingly, 
only  the  8n  and  8n-4  cases  are  of  practical  imiportance.  These  two 
cases  are  shown  in  Figs.  ?  and  h  for  n=2. 

The  parameters  e,  s,  and  b  will  be  computed  below  for  the  8n  and 
8b-4  sided  polygon  cases.  The  subscript  p  will  be  used  In  the  8b-^ 
case  equations  to  prevent  confusion. 


fricliwpmMyt 


The  sides  of  the  regular  polygons  are  defined  by  the  angles  of 
minimum  | w | , 

*ni  »  1  S  m  s  n  (8n  case)  ,  (3-1) 

♦q  "  ’  1  ^  q  ^  P  (Sp-^^  case)  ,  (5-2) 

the  vertices, 

^  )  0  ^  m  ^  n  (8n  case)  ,  (3-3) 

’•'vq  "  ^^i-2^”  '  1  ^  q  ^  P  (8p-^  case)  ,  (3-4) 

and  the  fact  that  the  vertices  are  on  the  unit  circle  in  the  8n  case, 

ik)  ■  <5-5) 

and  the  sides  are  tangent  to  the  \init  circle  in  the  ^-4  case, 

aj  +  b^  =  1  .  (5-6) 

Using  Eqs.  (3-5)  an<i  (5-6)  with  b/a  =  tant,  where  y  is  given  by 
Eq.  (3-1)  and  (3-2),  yields 

^m  =  [*  ^  y 

=  r  cos(e-yjij)  sec^^^  ,  (3-7) 

l4 


and 


f ^  =  X  cos(\|r^)  +  y  sin(y^) 


=  r  coB(e-if^) 


(3-8) 


For  the  8n  case,  the  peak  error  is  at  the 


f(t  )  > 

-n  =  -T-"  ■  ^  ^  ^ 


(3-9) 


■while  for  the  8p-4  case,  the  peak  error  is  at  the  vertices. 


e  =  1  - 
P  r 


H. 


=  1  -  cos 


(3-10) 


The  multiplicative  bias  is  found  by  integrating  r/r  with  respect 
to  9  over  one  side  of  the  polygon  and  dividing  by  the  angle  between 
the  -vertices,  since  9  is  \aiiformly  distributed  over  the  sector  corres¬ 
ponding  to  each  side  of  the  w  plane  polygon. 


Thus, 


n  ^  ^  y  ^  »  ^5-11) 

^  o 


* 

BE 


and 


P  * 


L 


*/(%>-M  SlP 


cos9(t9 


X 

apr 


(3-12) 


using  the  obvious  expedient  of  integrating  over  half  a  aide. 
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(3-16) 


^.0 


Implementation 


IKe  Molest  ..y  to  E,.  (3.7,  „  j,,  (3.3, 

*1  ■  ‘“(5;)  (a»  '^«e)  (t.l) 


or 


=  X  (^>.4  case) 

Then,  fot  1  .  .  .  „.i  „  1  ^  ^ 

*m+l  *  *'ni  “  ^  +  y  ®®®(^)) 


or 


Vl  ”  ^q  =  2  sln/^ 


-  ®‘"'^)  *  y 


can  be  used  to  fonn  9^,  ?  ...,  as  follows. 


Conqjute 


A  * 
m 


y  •  X  ti 


-(g) 


(^-2) 


(M) 


or 


(^-5) 


(4-6) 
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^"*•*•1  «•  Ihrt 


If  A  ^  0,  the  computation  is  completed;  otherwise,  use 


r  -  r  =  2  sin 
q+L  q 


(4-8) 


to  update  f  or  f  .  If  m^n  or  q=p,  the  con?)utation  is  finished; 
m  q 

otherwise,  repeat  the  confutation  of  A  and  subsequent  sign  check. 


The  quantities 


Cj^(tn)  =  tan^^l  ,  1  5  m  ^  n-1 


>)  =  2  tan(Jj)  cos(g) 


1  s  m  s  n-1 


(4-10 


c^(q)  =  tan 


1  s  q  S  p-1 


(4-11 


C2(q)  =  2  sin^^^j^  cosi 


,  1  S  q  S  p-l  ,  (4-12 


along  with  taii(«/8n),  e  ,  e  ,  b„,  b  ,  s  ,  and  s  ,  are  given  in  the 
'  '  n  p  n  p'  n  p 

appendix,  for  all  n  and  p  from  1  to  12,  in  both  decimal  and  octal. 

The  nniltiplications  in  Eqe.  (4-l),  (4-5),  (4-6),  (4-7),  and  (4-8) 
can  be  accoinplished  within  the  purposes  of  the  algorithm  by  using  only 


the  top  L  bits  of  tan(mn/4n),  Cj^  and  where  L  is  the  niariber  of 
zero  leading  bits  In  b  or  S,  depending  on  whether  neurrowband  flutter 
or  slgnal-to-nolse  Is  the  limiting  factor.  In  the  CDC  3200,  if 
(A)=-y  and  (Q)a+x,  -A  can  be  computed  faster  by  several  Ehift>add 
(SI^,  AQ/V)  Instructions  than  by  a  single  multiply-add  (M(JA,  ADA)  execu¬ 
tion  if  the  number  of  nonzero  bits  among  the  most  significant  L  In 
the  multiplier  is  less  than  six.  An  AZJ,QE  Instruction  can  then  be 
used  to  retrieve  r  and  exit  If  A  §  0.  ;nie  p=2  (regular  12  sided 
polygon)  case  Is  progranned  as  an  exazQile  In  the  appendix;  entering 
with  X  and  y  in  (a)  and  (Q),  the  execution  times  eu^e  17  usee  minimum, 

8  S 

2k  usee  average,  and  3^  Msec  meudmum,  including  call  euid  return  with 
the  result  in  (A).  This  routine  is  on  disc  file  SPAM. 
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5.0  Appendix 


5.1  Exaagle 

This  is  a  listing  of  a  COMPASS  routine  for  the  12  sided  polygon 
envelope  estitnate  shown  in  Fig.  3*  The  peak  error  is  -K),  -0.054^^^, 
the  multiplicative  bias  is  O.989,  and  the  rms  relative  error  is  0.01^^^ 
so  that  narrowband  flutter  is  O.15  dB  peak-to-peak  and  the  maximum 
signal -to-noise  is  4o  dB.  The  execution  times  are  I7  psec  minimm, 

24  psec  average^  and  3^  Psec  maxlmtm,  including  call  and  return. 

Entry  is  accomplished  by  a  BTJ  XY.SQRT  instruction  with  the  fixed 
point  quadrature  sauries  in  (a)  and  (Q),  and  return  is  to  P+1  with  the 
envelope  estimate  in  (a). 


rnoAifPlillMk 
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5.2  Tables 


Following  are  tables  of  constants  to  be  used  in  choosing  and 
Implementing  an  envelope  estimate  based  on  a  Uk  sided  polygon.  For 
any  k,  the  peak  relative  error  is  -fO, -EP  for  odd  k  or  +EP, -0  for  even  k, 
the  rms  relative  error  is  20  log^^SP,  and  the  multiplicative  bias  is 
BP.  If  L  is  the  number  of  leading  zero  bits  in  BP  or  SP,  implementation 
is  as  follows. 

Find  X  =  absolute  value  of  largest  sanple,  y  =  absolute  value  of 
smallest  sample.  Define  an  initial  r  by 

r  =  X  ,  k  odd  ;  r  =  x  +  C^y  ,  k  even  .  (5-1) 

Use  shift-adds  for  the  multiplication,  using  only  the  nonzero  bits  in 
the  L  bits  of  following  the  binary  point.  Then,  compute 

-A  =  -y  +  x*Cl(l)  ,  (5-2) 

using  shift -adds  with  the  nonzero  bits  among  the  L  most  significant 
in  Cl(l)  for  the  multiplication.  If  -A  *  0,  exit  with  r  in  (A). 
Otherwise,  vpdate  r  using 

?  =  r  +  A*C2(1)  ,  (5-3) 

using  shift-adds  with  the  nonzero  bits  among  the  L  most  significant 
in  C2(l)  for  the  multiplication,  as  before.  If  no  Cl's  or  C2'b  reosdn, 
exit  with  9  in  (A).  Otherwise,  repeat  the  qperatioos  and  checks 
beginning  with  Eq.  (5-2)  using  Cl(2)  and  C2(2),  Cl(5)  and  C2(5),  etc., 
until  either  ^  k  0  or  there  are  no  laore  Cl* a  or  C2'b. 
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